New limit insertion on pseudoscalar-photon mixing from WMAP Observations 
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The pseudoscalar-photon mixing in presence of large scale magnetic field induces polarization in 
light from distant cosmological sources. We study the effect of these pseudoscalars or axion like 
particles (ALPs) on Cosmic Microwave Background Radiation (CMBR) and constrain the product 
of mixing strength g^ times background magnetic field B. The background magnetic field has 
been assumed to be primordial and we assume large scale correlations with the correlation length of 
IMpc. We use WMAP seven year foreground reduced polarization and temperature data to constrain 
pseudoscalar-photon mixing parameter. We look for different mass limits of the pseudoscalars and 
find g^B < 1.6 x KT 13 Ge V~ x nG with ALPs of mass 10~ 10 eV and g<t> B < 3.4 x IQ^GeV^nG for 
ultra light ALPs of mass 10~ 15 eV. 



I. INTRODUCTION 

The pseudoscalar-photon mixing and its effect on the 
distant cosmological sources has been studied in lit- 
erature [H-l2lj|. These hypothetical axion like parti- 
cles(ALPs), arise naturally as pseudo-Goldstone bosons 
in theories with spontaneously broken global symmetries 
j22Tj30j . ALPs have an interaction vertex with two pho- 
tons and hence in an external magnetic field ALPs can 
convert into a photon and vice versa (3lTl4^ . Although 
this effect is very small, it becomes significant at cosmo- 
logical scales and leads to many interesting signatures 
on electromagnetic radiation. This pseudoscalar-photon 
mixing phenomena causes changes in intensity as well 
as polarization in radiation from distant sources (ll-l2l|. 
The contribution of this effect has been investigated for 
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sources. Various experiment are looking for these pseu- 
doscalars and providing limit on the coupling constant 
50 and their masses m^H, [IH, l5lTi70| . 

In the present paper we study the effect of 
pseudoscalar-photon mixing on CMBR multipolc 
anisotropy. We show, using WMAP observations that 
this leads to a new constraint on the product of mag- 
netic field B and the pseudoscalar-photon coupling g^. 
We consider the background as the large number of 
correlated magnetic field domains and do a complete 
3-D-simulation [l6[ to calculate the stokes parameters 
for CMBR. The origin of background magnetic field 
is considered as primordial [7ll - [75| and we assume a 
smooth variation of the magnetic field over the scale of 
1 Mpc. However, as we do the simulation over a very 
large distance (redshift 1000), we choose the domain size 
around 16Mpc and the strength of magnetic field in each 
domain is order of nG @, S HH • We compare our result 
with the WMAP 7- year data and constrain the coupling 
parameter g^ times B. 



The paper is organized as follows. In Sec|TI]we briefly 
review the pseudoscalar-photon mixing in presence of 
plasma and uniform magnetic field in a flat expanding 
universe. In Sec lIIII we model the background magnetic 
field, which is correlated in real space and discuss the 
numerical method for generating the ZD magnetic field. 
In Sec lIVI we present our simulation result and compare 
with the WMAP observations. Finally, in Sec|V]we con- 
clude and compare our results with available literatures. 



II. PSEUDOSCALAR-PHOTON MIXING 

A. Basic Formulation 

In this section we briefly describe the propagation of 
electromagnetic waves. We choose a fixed coordinate sys- 
tem such that the z-axis lies along the direction of prop- 
agation. We define A. = ^ E ' where E is the usual elec- 
tric field vector and u is the radiation frequency. Only 
the component of A parallel to B transverse mixes with 
the pseudoscalar field 4>. We replace <fi by ^ and the 
mixing of A\\ to \ can be written as, 



— M 



= 0, 



(1) 



where M is the 'mixing matrix' as, 
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Here m^, is the pseudoscalar mass and (a 2 B±) is the 
transverse component of magnetic field, the factor a 2 is 
scaling the magnetic field in expanding universe model. 

We follow the procedure described in Ref. [4(| EH for 
solving Eq[TJ 



2 



B. Propagation and Polarization as, 



In this section we briefly review the propagation of the 
mixed A and x field in presence of an external magnetic 
field and calculate Stoke's parameters to quantify polar- 
ization. A detailed description is given in Ref.fiol l43j . 
We start with some initial densities of A and \ mixed 
fields and calculate the same after a distance z propaga- 
tion. An appropriate and general representation is given 



P(0) 

and 

P(z) 



where 7 = 1 — e lAz , while A = — A^ and A^ , 
A^ are defined in terms of cj and eigenvalues, of the 
matrix M, 



Aa = ^u 2 - \i\ - uj, = ui 2 — p?_ lo. (6) 

We assume that the CMBR is unpolarized at z = 1000 
and also set initial density for pseudoscalars to be zero. 
Initially the plasma density is supposed to be very low 
(n e = 3.24 x 10~ 10 a- 3 cm- 3 ) [7^ (z re 1000) as the uni- 
verse has gone through the recombination (z re 1100) era 
and almost all the electrons and protons have been com- 
bined to form neutral hydrogen. At redshift z re 6 [T7| . 
the star formation starts and the universe becomes ion- 
ized, so for redshift z < 6 we assume plasma density 
n e = 10~ 8 a~ 3 cm~ 3 . We next propagate the CMBR 
through these domains and compute p{z) in each domain. 
Finally we obtain p(z) and calculate the Stoke's param- 
eters as, 

I(z) =< A\\{z)A\{z) > + < A±(z)Al(z) > (7) 
Q(z) =< A\\{z)A\{z) > - < A±(z)Al(z) > (8) 
U(z) =< A l{ (z)A* ± (z) > + < A±{z)A\{z) > (9) 
V(z) = i(- < A\\(z)A* ± (z) > + < A±{z)A\{z) >) (10) 

We compute multipole anisotropy in E and B modes and 
constrain g^B by demanding consistency with CMBR 
observations. 



III. BACKGROUND MAGNETIC FIELD 

It is reasonable to assume the orig in of background 
magnetic field as primordial [7ll - [75l l78j| . A two-point cor- 



p(z) = P(z)p(0)P(z)-\ (3) 

where p(0) is the general density matrix for the fields 
and P(z) is a non trivial unitary matrix, describing the 
solution for the mixed fields. The representation for p(0) 
and P{z) for our coordinate system is as follows, 



(4) 



(5) 



I 

relation function for a homogeneous and isotropic mag- 
netic field is given as, 

(b i (k)b*(q))=5 k , q P ij (k)M(k) (11) 

where bj(k) is the j th component of the magnetic field in 
wave vector space. The real space magnetic field Bj(r) 
can be written as a Fourier transform of fcj(fcJE 

Here 

Pij{k) = — ^j^J is the projection operator and 
function M(k) is given as, 

M(k)=Ak nB ; (14) 

where tib is power spectral index and the constant A 
is a normalization. The numeric value for A is such as 
X)< < Bi(r)Bi(r) >= B 2 , B is the strength of the 
magnetic field, often assumed to be InG |73l. |79| on a 
comoving scale of IMpc. We assume ub = —2.37, which 
is derived from the best fit of matter and CMBR power 
spectrum [79|]. Our results depend on jib and a different 
value will change the imposed limits. We do not impose 
any cutoff on correlations in real space and simply choose 
r max larger than our system. In other words the lower 
limit on wave vectors k m i n = rzS tends to zero. 



The real and wave vector space field transformation are as fol- 
lows, 

^w^Ey*"^ ( 12 ) 



<^H (0)^(0) > < A\\(P)A*j_(P) > <-4||(0)x*(0)> 
<A±(0)AI(0)> < A±(0)Al(0) > <A±(0) X *{0)> 
<X(0MI(0)> <x(0)A* ± (>)> <x(0)x*(0)> 



1 — ^jsin 2 6 ^jcosOsinO 

e i(cu+A A )z I Q e _ i[w+A<1 -( tl , 2 -a,= ) 1 /2] x Q 

jcos8sin8 1 — jcos 2 9 
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We split the space in 1024 x 1024 x 1024 equal volume 
domains and generate the 3D k-space magnetic field in 
each domain using the spectral distribution as in EqfTTI 
We use polar coordinate (k, 9, <J)) in wave vector space. In 
K-space the domains are uncorrelated and for any wave 
vector k, bk = and bg and b^ are uncorrelated. Hence, 
we can generate the bg and b^ independently for each 
domain using a smooth Gaussian distribution|l6l. l46j . 



f(b e {k),bt{k))=N exp 



2M(k) 



,(15) 



Here N is a normalization factor. We use this to gener- 
ate full 3D k-space magnetic field for each domain and 
do a Fourier transformation to get the three Cartesian 
components of the magnetic field in real space. 



IV. SIMULATION AND RESULT 




10 z 10" 1 
scale factor(1/1+z) 



FIG. 1. The mixing angle [&) in flat expanding universe, 
m = 10~ 10 eF, g4 ,B = 1.6 x W^ 1 ' 3 GeV^nG. 



We propagate CMBR from redshift 1000 and do the 
simulation for the pseudoscalar-photon mixing in the cor- 
related magnetic field background. We perform our com- 
putation on a 3D grid of 1024 x 1024 x 1024. Since the 
total linear distance for redshift 1000 correspond to a very 
large comoving distance 8104Mpc(matter dominated uni- 
verse), we set one domain size to be 16Mpc however we 
still keep the correlation length fixed to IMpc. This do- 
main size is very large as compared to the oscillation 
length (0.4Kpc) for the CMBR and it has been observed 
that a small domain size leads to a strong mixing effect 
if it is larger than the oscillation length. However we do 
not expect a very significant effect due to the domain size. 
We have perform simulation with different domain size 
in Ref. (l6l | and observed very small variation statistically. 
We use HEALPis0 to generate angular positions of the 
sources with resolution parameter 256(pixel size of 13.7 
arcmin). Next, we propagate all these sources through 
each domain and determine the Stoke's parameters. 

We compare our simulations to seven year foreground 
reduced WMAP CMBR observations. We use W-band 
since it contains the least foreground contamination. We 
demand that the pseudoscalar coupling introduces the 
temperature fluctuation less that 10~ J . Furthermore we 
impose the constrain from E and B mode and demand 
that the mixing mechanism generate E and B mode 
within the limit of observation. In our analysis we 
always have the coupling constant g^ multiplied with 
the background magnetic field B and so we are able 
to put a limit only on the term g^B. We simulate the 
CMBR polarization with the following parameters:- 

(1) B T = InG 

(2) plasma density n e — 10~ s a~ 3 cm~ 3 for z < 6 and 
n e = 3.64 x 1Q- 10 a- 3 cm~ 3 for z > 6 [76| 



2 http://hcalpix.jpl.nasa.gov/ 
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scale factor(1/1+z) 



FIG. 2. The mixing angle(#) in flat expanding universe, 
m = W~ 15 eV, g<t,B = 3.4 x 10" 15 GeF _1 nG.Here plasma 
mass term(Ljp / a) is dominating over pseudoscalar mass term 

K« 2 ). 



(3) CMBR frequency for W-band = 90GHz 

(4) pseudoscalar mass =10~ 10 eV and 10~ 15 eV 

It turns out that most of the polarization and 
anisotropy is generated at high redshift. At low redshift 
the mixing is negligible as shown in Fig. (|ll2[ ). We sim- 
ulate the mixing effect for = 10 _10 eV and = 
10- 15 eV. 

We have clearly seen that the pseudoscalar mass is a 
very significant parameter and the predicted CMBR po- 
larization increases rapidly with decreasing pseudoscalar 
mass unless it touches the plasma mass limit. 

CMBR studies has a convention to use spherical har- 
monics and drive the multipoles for anisotropy studies. 
We follow the same and produce a full sky map of Stoke's 
parameters with 13.7 arcmin resolution for W-band fre- 
quency and calculate E and B modes. We degrade the 
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E mode multipole, = 10 10 eV, g^B — 1.6 x 

W~ 13 GeV~ 1 nG FIG. 5. The simulated(red) and WMAP observed(blue) 

E mode multipole, = 10~ 15 eV, g^B = 3.4 x 

lO' 16 GeV^nG. 




FIG. 4. The simulated(red) and WMAP observed(blue) 
B mode multipole, = 10~ eV, g^B = 1.6 x 

\Q- 13 GeV- 1 nG. 




FIG. 6. The simulated(red) and WMAP observed(blue) 
B mode multipole, = W~ 15 eV, g^B = 3.4 x 

W^GeV^nG. 



WMAP seven year foreground reduced I,Q,U sky maps 
to match with our simulation resolution and deduce the 
E and B modes for W-Band frequency. We choose an 
appropriate value for the factor g^B requiring the E and 
B mode multipoles within the observation. This leads 
to a limit on the factor g^B. We have shown our re- 
sults in Fig. dSHJl for m = lCT 10 e^ and in Fig.flSD for 
= lCP 15 eT^. If we set = 10~ 10 eV the limit is fixed 
to be g^B < 1.6 x KT^GeV^nG. Alternatively if we 
choose ultra light pscudoscalar of mass < 10~ 15 eV, 
then we fixed g^B < 3.4 x 10~ 15 GeV~ 1 nG. We present 
the temperature multipole anisotropy in Fig. ()7I5|) . the 
upper curve(blue) is the WMAP observation and the 
lower(red) one is our simulation. We find that with our 
limits, the simulated temperature anisotropy is below 
the WMAP data. Our constrain on pseudoscalar-photon 
mixing are deduced from the E, B modes of CMBR ob- 
servations. 



V. DISCUSSION 

Wc have done full 3D simulation of pseudo -scalar pho- 
ton mixing for CMBR at a very high resolution of 13.7 
arcmin over the full sky. A comparison with WMAP ob- 
servation results as a new and more stringent limit on 
the factor g^B. It depend on the pseudoscalar mass and 
we simulate the limit on the factor g^B for two different 
masses of pseudoscalars. 

Recently [12], a bound on factor g^B as g^B < 
10~ 11 GeV~ 1 nG has been derived from ultraviolet photon 
polarizations emerging from active galactic nuclei. Here 
the derived limit corresponds to ultra light ALPs(m^, < 
lCT 15 eV^). In Ref.flU l80j the limits on g^B has been 
studied through CMBR spectral distortion, giving g^B < 
10~ 13 ~ 10~ u GeV _1 nG for ALPs masses between 
lO -15 ^ and 10~ 4 eV. The pseudoscalar-photon mix- 
ing may also contribute to the dimming of Type la 
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FIG. 7. The simulated(red) and WMAP observed(blue) 
T mode multipole, = W~ 10 eV, g^B = 1.6 x 

10- 13 GeF _1 nG. 




FIG. 8. The simulated(red) and WMAP observed(blue) 
T mode multipole, — 10 _15 eV, g^B = 3.4 x 

10" 15 GeF _1 nG. 



supernovaejl, 0, [Slj] . The phenomenon fixes g^B to 
10- 11 GeF _1 nG for a axion of mass 1(T 16 ey[| 

We may possibly constrain g^ form our bound on g^B, 
however the very uncertain large scale background mag- 
netic field B leads some uncertainties. Assuming the 
background magnetic field as InG over a comoving scale 
of ps lOMpc, our results bound g$ < 1.6 x IQ-^GeV' 1 
and g<j> < 3.4 x lO^GeF" 1 for the ALPs of \Q- w eV 
and 10 _15 ey respectively. 

Our limits can be compared with the direct experimen- 
tal limits from SN1987A , which is 50 < lO^GeV^ 
and ,g < 3 x \Q- 12 GeV^ for very light ALPs (< 
10~ 9 eV). We also recall the results from CASTA|g!|67j|, 
,g < 8.8 x KT n GeV for the ALPs of 0.02eV, which of 
course not for the ultralight ALPs and do not have a 
direct comparison with our results. 

We conclude that the CMBR multipole anisotropy im- 
poses a stringent constraint on the pseudoscalar-photon 
coupling. Wc have obtained the lowest value of g^B as 
compared to available literatures. 
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